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Abstract 

The Kumaraswamy Inverse Weibull distribution has the ability to 
model failure rates that have uninrodal shapes and are quite common in 
reliability and biological studies. The three-parameter Kumaraswamy 
Inverse Weibull distribution with decreasing and unimodal failure rate 
is introduced. We provide a comprehensive treatment of the math¬ 
ematical properties of the Kumaraswany Inverse Weibull distribution 
and derive expressions for its moment generating function and the r-th 
generalized moment. Some properties of the model with some graphs of 
density and hazard function are discussed. We also discuss a Bayesian 
approach for this distribution and an application was made for a real 
data set. 

Keywords: Kumaraswamy distribution, Weibull distribution, Survival, 
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1 Introduction 


In this paper we propose a new probability distribution to handle the prob¬ 
lem of survival data. Motivated by research developed in recent years, we 
introduce the Kumaraswany Inverse Weibull distribution that includes sev¬ 
eral well known distributions used in survival analysis. 

Recently, many authors have proposed new classes of distributions, which 
are modifications of distribution functions which provide hazard ratios con¬ 
templating various shapes. We can cite for example the Weibull exponential 


flMudholkar et al 


19951 ) , which has a l so th e h azard rate func t ion w ith a 


unimodal form, (see also 


Xie and Lai (19 96 1)) 


Carrasco et al. 


(120081 ) pro¬ 


posed a four-parameter dist ribution denoted generalized modified Weibull 
(GMW) distribution, 


Gusrnao et al. 


(201 lj) introduced and studied the tri- 


parametric inverse Weibull generalized distribution t. 
rate with unimodal, increasing and decreasing forms. 


iat possesses 

failure 

Pescim et al. 

(2010) 


proposed a distribution with four parameters, called beta generalized half 
normal distribution. 

Underexplored in the l iterature and ra r ely us ed by statisticians, the Ku- 


maraswamy distribution ( Kumaraswamvl . 


19801 ) has a domain in the real 


interval (0,1). This property turns the Kumaraswamy distribution a natu¬ 
ral candidate to combine with other distributions to produce a more general 
one. Its cumulative distribution function (cdf) is given by, 

F(x;a,b) = 1 - [1 - z“]\ (1) 

and its probability density function (pdf) is given by, 

/ (x; a, b) = abx a ~ l [1 — x a ] b ~ 1 , 0 < x < 1, 
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where a > 0 and b > 0. This density can be unimodal, increasing, decreasing 
or constant. 


Recently, 


C ordeiro and de Castrol (201 lj) proposed to use the Kumaraswamy 


to generalize other distributions. Considering that a random variable X has 
distribution G, they suggest to apply the Kumarasawamy distribution to 
G{x). Note that, since 0 < G(x ) < 1 for any distribution function G, then 
evaluating equation (JTJ) at G(x) we obtain, 


F G (x-,a,b) = l-[l-G(x) a ] b . 


( 2 ) 


where Fq is the cdf of the generalized Kumaraswamy-G distribution. Based 
on these ideas, we consider the Inverse Weibull Distribution as a candi¬ 
date for G, using Equation Q. Then, performing some adjustments and 
mathematical manipulations, we obtain the Kumaraswamy inverse Weibull 
(Kum-IW) distribution. 

The rest of the paper is organized as follows. In Section^ we develop the 
Kum-IW distribution. Section [3] is devoted to describe basic properties of 
the distribution. Inference procedures via maximum likelihood and Bayesian 
approaches are presented in Section 01 Section [5] is devoted to analyze a real 
data set and in Section E we present some conclusions of this work. 


2 Kumaraswamy inverse Weibull distribution 


Let T a random variable with inverse Weibull distribution. Then its cdf can 
be written as, 


G(t ) = exp 



t > 0, 


( 3 ) 
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where a> 0, /5 > 0, and its pdf is given by, 


g(t ) = j3a.Pt 


exp 


a\P 
t 


Inserting the G function of Equation (|3|) in Equation (J2J) it follows that, 


Fa(t; a,b, a, (3) = 1 — < 1 — exp 


a\P' 
~ a 1 1 


( 4 ) 


We note that the parameters a and a in Q are not identifiable and we adopt 
the reparameterization c = aa 1 ^ so that the Kum-IW cdf is rewritten as, 

0-1 N b 


Fait ; b, c, f3) = 1 — < 1 — exp 


c\P 
t 


(5) 


where b > 0 and c > 0 are the shape and scale paremeters respectively. 
Accordingly, the Kum-IW pdf is now given by, 

o-| A 6-1 


f G (t\b,c,p) = pbcPt-V+V 


exp 


c\P 
t 


1 — exp 


c\P 
t 


( 6 ) 


It can be easily seen that when b = 1 we obtain the pdf of the Inverse 
Weibull (IW) distrbution given by, 


g{t) = j3aPr^ +l) 


exp 


-■!)' 


Finally, the corresponding survival and hazard functions are respectively 
given by, 

PbcPt-(P + V ex p [- (f) / 


S G (t]b,c,/3) = < 1-exp 


c\0 
t 


and h G {t\ b , c, j3) = 


1 — exp Y~ (f)' 

while the quantile function, Q(u), of the Kum-IW distribution is given by, 

Q(u) = F~ l (u] b, c, f3) = c(-log(l - (1 - 
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2.1 Some special classes of the Kum-IW 

The following well known and new distributions are special sub-classes of 
the Kum-IW distribution. 


• Kumaraswamy inverse Rayleigh distribution (Kum-IR) 


If /3 = 2, the Kum-IW distribution reduces to the Kumaraswamy 
inverse Rayleigh distribution (Kum-IR). Then, with j3 = 2 the density 
function of Kum-IW is expressed by: 


F G (t\b,c ) = 1 - < 1 - exp 


c\2 

t 


where 6 > 0 is the shape parameter, and c > 0 is the scale parameter. 
Hence, the KiR distribution has two parameters, and its pdf is given 
by 


fc(t- b, c ) = 2bc 2 t 3 exp 



1 1 — exp 




The corresponding survival and hazard functions are given respectively 
by, 


S G (t;b,c ) = <M - exp 


and h G {t-,b, c) = 


2bc 2 t 3 exp —(f)" 


1 — exp 


-(f)' 


Inverse Rayleigh distribution (IR) 

If (3 = 2 and 6=1, the Kum-IW distribution reduces to the inverse 
Rayleigh distribution (Kum-IR). Then, with (3 = 2 and 6 = 1 the 
density function of Kum-IW is expressed by: 

2 " 


G(t) = exp 


t > 0, 
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where a > 0, and its pdf is 


g(t) = 2 a 2 t 3 exp 




Kumaraswamy inverse Exponential distribution (Kum-IE) 

If /3 = 1, the Kum-IW distribution reduces to the Kumaraswamy 
inverse Exponential distribution (Kum-IE). Then, with (5=1 the 
density function of Kum-IW is expressed by: 

F c (t:b,c) = 1 - jl - exp - (^) } , 

where 6 > 0 is the shape parameter, and c > 0 is the scale parameter. 
Hence, the KIE distribution has two parameters, and its pdf is given 
by 

—2 


f G (t]b,c) = bet 2 exp 


©K—KDr 1 ' 

The corresponding survival and hazard functions are respectively 
S c (t:b,c) = jl - exp ^ 

Inverse Exponential distribution (IE) 


c\iI b bet 2 exp [— (§)1 

t) } and h G (t;b,c) = - - f (c ^ 

tn) ! — exp [— (|)J 


If /3 = 1 and 6=1, the Kum-IW distribution reduces to the Inverse 
Exponential distribution (IE). Then, with (3 = 1 and 6 = 1 the density 
function of Kum-IW is expressed by: 


G(t ) = exp 
where A > 0 and its pdf is 

g(t ) = A t~ 2 exp 


t > 0, 
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(a) (b) 

Figure 1: Kumaraswamy inverse Weibull density functions. 

3 Basic Properties 


In this section we describe in detail some properties like expansions, mo¬ 
ments, mean deviations, Bonferroni and Lorenz curves, order statistics and 
entropies which might be useful in any application of the distribution. 


3.1 Expansions for the distribution and density functions 

We now give simple expansions for the cdf of the Kumaraswamy Inverse 
Weibull distribution. If |x| < 1 and ip > 0 is a non-integer real number, we 
have 

OO 

(1-*)* = £(-1)‘«>!)z‘. (7) 

7=0 

If ip is a positive integer, the series stops at i = ip. Using expansion in 
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(a) (b) 

Figure 2: Kumaraswamy inverse Weibull survival functions. 


Equation (0 it follows that, 




= RhrPt-v+v y 


i =0 


(*■) r (6-<) exp 


-G)Vd 


( 8 ) 


and 


F(t;b,c,f» = l-gizEIMiexp 


£j(i!)r(&+l-i) 


. c\P 

-"i 


Because the integrals involved in the computation of moments, Bonferroni 
and Lorenz curves, reliability, Shannon and Renyi entropies and other in¬ 
ferential results do not have analytical solutions, these expansions are nec¬ 
essary. 


3.2 A general formula for the moments of the Kum-IW 

We hardly need to emphasize the need and importance of the moments 
in any statistical analyses, especially in applied work. Some of the most 
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(a) 


(b) 


Figure 3: Kumaraswamy inverse Weibull hazard functions. 


important features and characteristics of a distribution can be studied using 
their moments (e.g. tendency, dispersion, skewness and kurtosis). If the 
random variable T follows the Kum-IW distribution, its k -th moment about 
zero is given by, 


(t k ) = f t k /3bc p t~^ +l) exp 


1 — exp 

[-Fill 

' ' Jo 







(-l)’T(b) 

—' T (b — r) r\ 

r =0 v 7 


(r + i)f -1 r ( 1 — ^ 


The moment generating function M(z ) of T for \z\ < 1 is 


,k =0 


^ T (b — r) r! 
r=0 x 7 



Hence, for \z\ < 1, the cumulative generating function of T is 

-z\,{-lYT{b) 


K (z) = log 


EE 

_k=0 r =0 


— be 1 
k\ r (6 — r)r! 


(r + iyp - 1 
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We note that it was necessary to use the expansions previously presented 
for the results of this section. 


3.3 Mean deviations 

The amount of scattering in a population may be measured by all the abso¬ 
lute values of the deviations from the mean or the median. If X is a random 
variable with Kum-IW distribution with mean fi = E[X] and median M, 
then the average deviation from the average and the average deviation to 
the median are defined respectively by, 


/*oo r oo 

di(X) = / \x — n\ f(x)dx and ^{X) = / \x — M\f(x)dx. 

Jo Jo 


Using the density of extended Kum-IW and given that, 


/»oo 00 

/ xf(x)dx = bc^2' 

J Li „ r\ 


(-l) r T(b) ( 1 \ t" 1 (p-l d 


^ T(6-r)r! ^1 + ^7 P ’ 11 P 

where 7(a, x) is the lower incomplete Gamma function, it follows that 

MX, - 


:(! + >-) , 


' T(b — r)r\ V 1 + r 
r =0 v ' v 


P V 


and 




T(b — r)r\ V 1 + r 

r=0 v ' v 


3.4 Bonferroni and Lorenz curves 


Bonferroni and Lorenz curves are widely applied not only in economics to 
study income and poverty, but also in other fields such as reliability, demog¬ 
raphy, insurance and medicine. 
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Let then p = E(X) e q = F~ 1 (p-,6), where is the inverse function 

of the cumulative function of a random variable X. Bonferroni and Lorenz 
curves are defined by, 


1 f q 1 fi 

B(p) = — / xf(x)dx and L(p ) = — / xf(x)da 
PH Jo H Jo 

Then, using the expanded density, 


rq 00 

/ xf(x)dx = bey 
r=0 


(-i) r r(6) / 1 \? -1 . 


—' T(6 — r)r\ \ 1 + r 


p-l c? 

~r , ¥ 


(1 + r) 


where T(a,a:) is the upper incomplete gamma function. Therefore, we have, 


the Bonferroni curve which is given by: 


be 




(—l) r T(6) ( 1 \P~ 1 (/3-l C P 


pp T(6 — r)r\ \l + r 


P > 


> —r (1 + r ) ; 


the Lorenz curve which is given by: 




a ^ T(b — r)r\ \ 1 + r 

r=0 v 


- 4-1 


r 


P -1 cP , 

_ 5 _ V (1 + r) )' 


3.5 Order statistics and Shannon entropy 


Let T\ :n < T 2 , n < • • • < T n:n be the order statistics obtained from the Kum- 
IW distribution. The random variable T r:n , for r = 1 , ,n, denotes the 
r-th order statistic in a sample of size n. The pdf of T r:n written as, 


= C r:n F(t ) r - 1 [1 - F{t)] n ~ r f(t), t > 0, 


where f(t) and F(t) are given by equations (O and © and C r:n = n\/[(r — 
l)!(n — r)!]. 
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The fc-th moment of the rth order statistic is 

OO 

/4$ = E(T k J = t k C r:n J [F^r 1 [1 - F(t)} n ~ r f(t)dt, 


for k = 1,2,..., and 1 < r < n. Hence, 


$l = C r:n [ 1 j n - r (l-j) n ~ 1 dj. 
J o 


and we obtain an expression for the moment given by, 


E(TL) = Cr.nhct Y.Y.1U 


(—l )*" 1 "- 7 T (bn — br + bj + b) 


-(z + l)* _1 T I 1 - 


j =0 i=0 


i!j!T (r — j) T (bn — br + b — i) 


The Shannon entropy of a random variable T is defined as a measure 
of the quantity of information. A certain message has more quantity of 
information the greater degree of uncertainty and is defined mathematically 
by E{— log[/(f)]}, where /(f) is the fdp of T. In particular, for a random 
variable T which follows the Kunr-IW distribution we have, 


E {- log[/(t)]} = - log (pbcP) +6(0+1) log(c) jr, r ( ft ) _r)H(r+ 1 ) 
i~(_iyr(&) 1 ™(-i yr(b), i 

^E r( 6-r)r! (r + 1) b ~ ^ + 1)1 + ° E f (h - rl rl (r + 1) ' 


i-l 


T (6 — r) r! 

r=0 x 7 


(—i) i+r r (b) r (6 + 1 ) .. . 2 . 

( 6 + EEjwr^rjrfj + i-j )^ 1 ) ’ J> 

r=0 i=0 v ' v 

where 7 = — log(j) exp {— j} dj is the approximate value of the Euler’s 
constant. 
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3.6 Renyi entropy 


The entropy of a random variable X with density function d2D measuring 
the uncertainty of the variation. The Renyi entropy is given by, 

I r ( p) = YZ~ p l0g {/ f ( X ) P dx } ' 


where p > 0 and p ^ 1. 

In information theory, Renyi entropy generalizes the Shannon entropy. 
This form of entropy is important especially in ecology and statistics, where 
it can be used as an index of diversity. In quantum information, it can be 
used as a measure of entanglement. If X is a random variable and follows 
the Kum-IW distribution, then the Renyi entropy is given by 


Ir (p) 



pp- l bc 1 ~ P Y^ 


i =0 


(-iyr(p b -p + i) 

(*!) r (pb - p +1 - i) 




i 

P 


4 Inference for Censored Data 


4.1 Maximum likelihood estimation 

Let Tj be a random variable with Kum-IW distribution with parameter 
vector 6 = ( b,c,/3 ). The data in survival analysis and reliability studies are 
generally censored. A very simple random censoring mechanism that is often 
realistic is one in which each individual i is assumed to have a lifetime T* 
and a censoring time Ci, where T, and Ci are independent random variables. 
Suppose the data set consists of n independent observations t, = min(T ? ;, Ci) 
for i = l,--- ,n. The distribution of Ci does not depend on any of the 
unknown parameters of Tj. Parametric inference for such data are usually 
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based on likelihood methods and their asymptotic theory. The censored 
log-likelihood 1(6) for the model parameters is 


£(6) = r log 


—c 




i€F 


i&F 


i£F 


1 — exp 


1 — exp 


( 10 ) 


6^ log 

ieC 

where F = [l,r] and C = [r + l,n]; still, C represents the censored data 
and F represents the failure data. 

The maximum likelihood estimate (MLE) 6 of d is obtained by solving 
the nonlinear likelihood equations Ub(6) = -j^p- = 0, U c (6) = = 0 

and Up(Q) = q„ ^ = 0. These equations cannot be solved analytically and 
statistical software can be used to solve the equations numerically. 

For interval estimation of b, c and /3, and tests of hypotheses on these 
parameters, we must obtain the 3x3 observed information matrix J(0) 
which is given by, 


( Jbb (o) Jbc (0) J b p (0) ^ 


J(9) = 


Jcb(0) Jcc(0) J c p(0) 

\ Jpb(0) Jpc(0) Jpp(9) J 

Under conditions met for parameters obeying the parametric space and 
not considering the limits of the same, the asymptotic distribution of 


V^(6-6) is N 3 ( 0 , T(6>)- 1 ), 


where 1(6) is the expected information matrix. This asymptotic behavior is 
valid if 1(6) is replaced by J(6), the observed information matrix evaluated 
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at 6. The asymptotic multivariate normal distribution N^(0, can be 

used to construct approximate confidence intervals and confidence regions 
for the individual parameters and for the hazard and survival functions. 
The asymptotic normality is also useful for testing goodness of fit of the 
three parameters the Kum-IW distribution and for comparing this distri¬ 
bution with some of its special submodels using one of the two well-known 
asymptotically equivalent test statistics - namely, the likelihood ratio (LR) 
statistic and the Wald and Rao score statistics. 


4.2 Bayesian approach 


Following the Bayesian paradigm, we need to complete the model speci¬ 
fication by specifying a prior distribution for the parameters. By Bayes 
Theorem, the posterior distribution is then proportional to the product of 
the likelihood function by the prior density. 

Subjetivism is the predominant philosophical foundation in Bayesian in¬ 
ference, although in practice noninf ormative prior densities (built on some 


formal rule) are frequently used (jKass et al 


1996). Since the parameters 


in the Kum-IW distribution are all positive quantities and due to the flexi¬ 
bility generated by the two-parameter Gamma distribution this is adopted 
as prior distribution. So, b ~ Gamma(m,w), c ~ Gamma(a, s) and j3 ~ 
Gamma(x, l). 

Assuming independence among the prior densities, the posterior density 
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is expressed by, 


h{0\t) OC h rn+n-l c a+nP-lpx+n-l 


exp 


—wb — sc — 1(3 — c^ ^2 


i =1 


if 

i =1 


-p-i 


1 — exp 


( 11 ) 


This joint density has no known analytical form but we can provide an 
approximate solution based on the complete conditional distributions of b, 
c and (3. These are given by the following expressions, 



n 

( n ( 

r / _ \ 

h(b c, (3, t ) oc b m+n 1 exp 

—wb — c p ^2 

\ IKH 1-exp 

i 

i 

1_ 


i =1 

[i=l l 


h{c | b. 


(3, t) oc c a+np 1 exp 


—sc — c 


i— 1 


n 

i= 1 




1 — exp 




n 

f n ( 


h((3 | b, c, t ) oc (3 x+n 1 exp 


) IIV j 1 -exp 

[- (2 ] 


i— 1 

[j=i l 


5 Application 


In this section, we present estimation results for parameters of the Kum- 
IW distribution under a Bayesian approach. The commercial production 
of cattle meat in Brazil, which usually comes from the cattle of the Nelore 
race, search to optimize the process trying to obtain a time for the cattle to 
reach the specific weight in the period of the birth until it weans. For a data 
set with 69 bulls of the Nelore race, was observed the time (in days) until 
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the animals reach the weight of 160kg relative to the period from birth until 
it weans. We compared the Kaplan-Meyer and Bayesian survival functions 
through two graphic methods. 


Using the expressi on o f 1(9 ) w 


Winbugs see 


Spiegelhalter et al. 


lich a routine was escribed in the software 


200 7) to esteem the values of the vector 


of parameters 6 = (6, c, (3) T of the Kum-IW distribution for the data of the 
cattle of the Nelore race. 


Table 1: Results of the Bayesian approach to Kum-IW distribution 


Parameter 

Mean 

SD 

2.5% 

Median 

97.5% 

b 

6.656 

11.18 

0.9829 

3.82 

29.76 

c 

177.5 

19.18 

154.9 

172.8 

227.5 

P 

8.231 

2.917 

3.84 

7.812 

15.01 


One of the methods that dispose for us to test if our model is well ad¬ 
justed to the data consists in the comparison of the function of survival of the 
parametric model proposed with the estimator of Kaplan-Meier. Another 
method consists of sketching the survival function of the model parametric 
versus the estimate of Kaplan-Meier for the survival function, if this curves 
is close of the straight line y = x we will have a good adjustment, figurc|5j 

6 Conclusions 

We worked a three parameter lifetime distribution called the Kumaraswamy 
Inverse Weibull (Kum-IW) distribution which extends Inverse Weibull dis¬ 
tribution proposed and widely used in the lifetime literature. The model 
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is much more flexible than the inverse Weibull. The Kum-IW distribution 
could have increasing, decreasing and unimodal hazard rates. We provide 
a mathematical overview of this distribution including the densities of the 
order statistics, Renyi entropy, Shannon entropy, Bonferroni and Lorenz 
curves and Mean deviations. Also, we derive an explicit algebraic formula 
for the r-th moment, expressions for the order statistics, and the maximum 
likelihood estimation for the censored data. The performance of the model 
was analized using real data sets where the Kum-IW distribution performed 
verywell and the estimation was given by Bayes method. 
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Figure 4: Trace plots of simulated parameter values of the Kumaraswamy 
inverse Weibull distribution. 
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(a) (b) 

Figure 5: (a) Comparison among survivals generated by the method 

Bayesian and Kaplan-Meier described in the graph S(t) versus time, (b) 
Comparison among survivals generated by the method Bayesian described 
in the graph S(t) versus Kaplan- Meier. 
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